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ABSTRACT 

The heating of the million degree, diffuse coronal plasma 
may be caused by a large number of events that are too 
small to be identified by present days observations. One 
explanation for these events could be the local interac- 
tion between magnetic flux systems that divide space into 
numerous flux regions. When such regions are indepen- 
dently advected by photospheric motions the expected 
outcome is the formation of enhanced current concentra- 
tion at specific locations in space. Due to magnetic re- 
sistivity, these currents dissipate and heat the plasma. In 
this paper, we investigate a simple model where two, ini- 
tially unconnected, flux systems are forced to interact in 
response to the imposed boundary driving by solving the 
non-ideal 3D MHD equations numerically. The recon- 
nection rate of the dynamical process is determined and 
compared with the corresponding rate for the potential 
evolution of the magnetic field. This shows that the dy- 
namic reconnection rate is about a factor of two smaller 
than the potential (perfect, instantaneous) rate for realis- 
tic solar driving velocities demonstrating that this three- 
dimensional magnetic reconnection process is fast. 

The energy input for a fixed advection distance is found 
to be independent of the driving velocity. The Joule dis- 
sipation associated with the reconnection process is also 
found to be basically dependent on the advection dis- 
tance rather than driving velocity. This implies that the 
timescale for the event determines the effect the heating 
has on the temperature increase. 

Finally, the numerical experiments indicate that the ob- 
servational structure of the reconnection site changes dra- 
matically depending on the phase of the evolution of the 
passage of the two flux sources. In the initial phase, 
where the sources become connected, the heating is con- 
fined to a compact region. For the disconnecting phase 
the energy gets distributed over a larger area due to the 
reconnected field line connectivity. 

Key words: Sun: Corona Heating, Magnetic Reconnec- 
tion. 

1. INTRODUCTION 

The coronal magnetic field is a manifestation of the dis- 
tribution of photospheric fluxes and their continuous mo- 
tion. By taking a magnetogram and representing each 
pixel value by a point source Close et al. (2004a) showed 
that each coherent flux patch connects to several neigh- 
bouring flux patches. This complex mapping of the mag- 
netic field gives rise to many separatrix surfaces; shells 



that divide space into different regions of connectivity. 
In general the distance, and also the direction, of travel 
of the various magnetic sources differ. The movement 
of any one source causes the orientation of the magnetic 
field at the interface of the separatrix surfaces enclosing 
this flux region to change. It is assumed that, in gen- 
eral, such change s lead to these surfaces attracting cur- 
rent dPriest et al. 2002) . If the current growth is signifi- 
cant magnetic reconnection is expected to occur releasing 
some of the free magnetic energy built up by the move- 
ment of the sources. 

In the photosphere the flux sources are continuously 
moved around by the convective motions and stresses like 
those indicated above must be taking place at many lo- 
cations above the photosphere. From detailed analysis 
of magnetograms, Hagenaar et al. (2003) showed that 
the recycle time of the photospheric flux is on the or- 
der of 8-19 hours. During this time flux sources emerge 
into the photosphere, where after they undergo merging, 
fragmentation and annihilation resulting in the loss of 
their identity. Flux sources are also advected resulting 
in a complex magnetic field. A direct consequence of all 
this activity is that the connectivity of the sources must 
change on a timescale that is much shorter than the tur n 
over time of the sources themselves dClose et al. 200 4b»). 
These processes all contribute locally to the heating of 
the coronal plasma, but their individual contributions may 
not be as prominent and easily observed as large flares. 
In the past, numerous models have been investigated for 
coronal heating. Most of these concern processes where 
flux emergence or annihilation releases the energy (e.g., 
Uchida & Sakurai, 1977; Forbes & Priest, 1984; Shibata, 
Nozawa & Matsumoto, 1992; Dreher, Birk & Neukirch, 
1997; Yokoyama & Shibata, 2001; Zhang & Low, 2002; 
Priest, Parnell & Martin, 1994; Parnell, Priest & Titov, 
1994; Birk, Dreher & Neukirch, 1996). In contrast, the 
situation where the flux is not emerged or annihilated, 
but instead is forced to interact with flux from a different 
source, have not been investigated in great detail. This 
process, which is highly likely to take place in the so- 
lar magnetic field, was first investigated using the min- 
imum current approach (Longcope 1998), where the in- 
teraction between independent flux sources in a constant 
background magnetic field is investigated. This work 
showed that as the flux systems start to interact the free 
energy in the field accumulates, in the form of a current 
sheet, at the separator line connecting two null points 
in the photospheric plane. Similarly, the same is found 
when the sources start to disconnect again. This scenario 
of independent flux systems first connecting and then dis- 
connecting was inv estigated using numerical experiments 
JGalsgaard et al. 2000 ). They showed that the initial pro- 



cess through which the two flux sources connect fol- 
lows the same pattern, namely that the two sources con- 
nect through separator reconnection. In the disconnection 
phase, however, the reconnection occurs through a differ- 
ent process - separatrix-surface reconnection. This stage 
of the evolution is clearly different from the minimum 
current approach. 

In this paper, we discuss the results of a series of numer- 
ical experiments that are conducted to investigate these 
basic processes of the flux interaction in more detail. In 
particular, the analysis concentrates on the reconnection 
process, determining a scaling relation between the nu- 
merical reconnection rate and the comparable rate for the 
potential evolution of the field structure, i.e., the perfect, 
instantaneous reconnection rate. Further, an analysis of 
the global energy input/output is investigated. Finally, a 
simple analysis of the magnetic structure of the newly 
reconnected field lines has been made, visualising the 
regions that would be exposed to large temperature in- 
creases as a result of the on going magnetic reconnection. 




Figure 1. An example of the initial conditions used in the 
experiments. The blue isosurfaces show the locations of 
the magnetic sources. The pink isosurfaces shows the lo- 
cation of two magnetic null points. The brownish field 
lines close to the top boundary shows the large scale 
magnetic field structure. The two other sets of field lines, 
red and green, indicate the location of the separatrix sur- 
faces dividing space into independent flux regions. The 
arrows at the bottom boundary indicate the profile and 
direction of the imposed driving velocity. 



2. MODEL SETUP 



The intention with the experiment is to investigate the 
simplest possible situation where two independent flux 
sources interact without loss (or gain) of flux from the 
sources. This is obtained by setting up a situation where 
two, oppositely-signed, flux sources are located a short 
distance apart on the bottom boundary of a Cartesian do- 
main. To remove the dipole structure the two sources 
would generate if on their own, a constant horizontal 
magnetic field is added. By changing the orientation of 
this overlying field component the field from the inde- 
pendent flux domains can be changed from running par- 
allel passed the other source, to parallel away from the 
other source. Introducing an advection velocity on the 
sources adds a preferred direction to the setup. The ad- 
vection velocity is applied such that the flux domains of 
the two sources are forced to interact with each other. 
As they interact strong current concentrations are gen- 
erated and eventually drive magnetic reconnection as the 
sources are advected passed one another In Fig. ^ the 
situation adopted in Galsgaard et al. (2000) is shown. 
Here, the isosurfaces show the locations of the magnetic 
flux concentrations. The field lines starting close to the 
bottom boundary indicate the separatrix surfaces divid- 
ing the space into three magnetically independent flux re- 
gions. In comparison, the large scale magnetic field is 
indicated by the field lines close to the top boundary of 
the domain. On the bottom boundary, the vectors show 
the structure of the imposed velocity. Notice, that the 
driving profile is such that it maintains the structure of 
the sources. 

The dynamical evolution of the system is investigated by 
solving the non-ideal 3D MHD equations numerically us- 
ing the code described by Nordlund & Galsgaard 1995. 



3. INTERACTION OF THE FLUX DOMAINS 

Experiments have been conducted with a number of dif- 
ferent orientations of the overlying magnetic field. The 
implications of this are, first, the distance the sources 
have to be advected before the two flux domains start to 
interact vary with the overlying field angle; and second, 
so does the velocity perpendicular to the field in the flux 
domains. At all orientations of the overlying field, the 
same general evolution takes place. As the two flux do- 
main are forces into each other a separator current sheet is 
formed connecting the two null points located on the bot- 
tom boundary. As the current magnitude reaches the nu- 
merical limit magnetic separator reconnection starts and 
new field lines are created connecting the two sources. 
This process continues until the two sources are well 
passed their point of closest approach and the two sources 
are well connected. Eventually current concentrations 
start forming over the separatrix surface that encloses the 
connected flux resulting in separatrix-surface reconnec- 
tion that disconnects this field. The two evolutionary 
phases are shown in Fig. |2] where the large isosurfaces 
represent the locations of the current sheets responsible 
for the reconnection and the field lines indicate the local 
and global connectivity of the magnetic field. 

4. MAGNETIC RECONNECTION RATES 

An important factor of magnetic reconnection in three- 
dimensions relates to the speed at which the process can 
take place. From two-dimensional magnetic reconnec- 
tion (see Priest & Forbes 2000 and references there in) it 
is know that the speed of magnetic reconnection depends 
critically on the boundary conditions of the investigated 
domain. Furthermore, most of the previous works have 
concentrated on steady-state situations. Here, the situa- 
tion is different, the system is driven and it reacts to the 



Figure 2. Left, the situation where separator reconnec- 
tion dominates the dynamical evolution. Right, the situa- 
tion at a later time where the field has started to discon- 
nect through separatrix-surface reconnection. The green 
isosurfaces indicate the locations of strong current, while 
the field lines indicate the structure of the magnetic field. 
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Figure 3. The field line connectivity at three differ- 
ent times. The colour coding represents different types 
of field line connectivity: blue indicates initially uncon- 
nected fiux, red represents connected fiux and green rep- 
resents the disconnected fiux. 



stresses created by the boundary motions and therefore 
finds it's own "boundary conditions" in the language of 
Priest & Forbes (1986) and in this sense finds its own 
steady state. 

To investigate the rate of change in connectivity a large 
number of field lines have been traced f rom the foot 
points of each source, (Parnell & Galsgaar^ 20041) . As 
the source structure is maintained in time, each field 
line can simply be associated with a fixed amount of 
flux. Therefore as the field lines change connectivity the 
amount of flux involved in each process can be counted. 
In Fig. |3l three snapshots showing the structure of the 
field line connectivity are shown. In these snapshots, 
blue represents the initially unconnected flux, red indi- 
cates the field lines that connect the two sources and the 
green area is associated with field lines that have been 
disconnected. These graphs show that the change in con- 
nectivity takes place in a non-trivial manner, providing 
a complex mapping of the field line connectivity. In the 
potential extrapolation of the coronal magnetic field by 
Close et al (2004a), it was shown that the field line map- 
ping of individual sources can connect to several neigh- 
bouring sources at the same time. From the mapping dia- 
grams shown in Fig.|3lit is seen that, even for this simple 
case, the connectivity map contains three different types 
of field line connectivity. This suggests that the mappings 
of Close at al (2004a) provide only a lower limit in terms 
of complexity and that the added effects of the dynamical 
evolution may increase the complexity of the field line 
connectivity of the coronal field even further. 

Another important effect of the dynamical evolution is 
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Figure 4. Left, the change in connected fiux as a function 
of time for the experiments with different angles of the 
overlying horizontal magnetic field. Right, the compara- 
ble graphs for the potential evolution of the experiments. 



the fact that multiple reconnection events take place at 
different spatial locations. This is clearly seen from the 
graphs showing the change in connectivity as a function 
of time where it is found that the reconnection processes 
connecting and disconnecting the flux systems, for a lim- 
ited time interval, take place simultaneous. In Fig.|4]an 
example of the time depended evolution of the connected 
flux is shown. The left panel shows the evolution in the 
situation where the magnetic field is assumed to always 
relax down to a potential magnetic field configuration. In 
the right panel, the comparable evolution is shown for 
the dynamical MHD development of the system. Since 
both connecting and disconnecting take place at the same 
time, to derive a simple reconnection rate for the system 
both processes must be taken into account. The process 
of tracking the field line connectivity provides a possibil- 
ity for doing exactly this. The amount of flux belonging 
to any of the three classes are therefore easily obtained 
as a function of time. Taking the time derivative of these 
and normalising the result with regard to the characteris- 
tic Alfven velocity of the domain and the total flux in a 
source, gives a dimensionless measure of the reconnec- 
tion speed. 

From Fig.l^it is seen that the fraction of connected flux 
between the two sources depends both on the angle of 
the overlying magnetic field and the driving speed and 
clearly on whether the evolution is potential or the dy- 
namical MHD solution. The time derivative of these 
curves gives a measure of the rate of change of the con- 
nectivity, which is equivalent to measuring a connect- 
ing reconnection rate. Similar graphs can be obtained 
for the disconnecting process and hence the disconnect- 
ing reconnection rate. By comparing these phases it is 
found that the connecting and disconnecting of the field 
in the potential case occurs at the same rate. In the dy- 
namical experiment, however, these processes are found 
to evolve at two different rates; with the peak connect- 
ing rate found to be about twice that of the peak discon- 
necting rate. This difference in reconnection rates is a 
consequence of the apparent stress of the magnetic field 
for the two process. In the connecting process the two 
independent flux systems are forced into each other as 
the system responds to the advection of the sources and 
the only way to release this stress is to generate a local 
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Figure 5. The peak total reconnection rate for the poten- 
tial evolution (dotted line) and the dynamical evolution 
(dashed line). The symbols represent the various individ- 
ual experiments while the lines are lines of best fit given 
by Eq. (0 and Eq. 0, respectively. 



current sheet through which magnetic reconnection can 
take place. In the following disconnecting process, re- 
connection takes place between open overlying magnetic 
field and the connected field which is advected under this. 
Here, the stresses between the two flux regions are less, as 
the overlying open field can simply be lifted up over the 
connected (dipole) field. The developing separatrix cur- 
rent sheets are consequences of the increasing tangential 
change in the field line orientation between these two flux 
systems and the pushing of them together by the magnetic 
tension and magnetic pressure forces in the two systems. 

Adding the rates of connecting and disconnecting gives 
the total reconnection rate at any one time. The result 
of this is seen in Fig. [Sj which shows both the potential 
(perfect, instantaneous) reconnection rate and the com- 
parable result for the dynamical evolution for all the con- 
ducted experiments. This graph doesn't provide detailed 
insight to the reconnection process, but shows how that 
the accumulated peak reconnection rate for the experi- 
ments depend, in a systematic way, on the driving veloc- 
ity perpendicular to the direction of the overlying mag- 
netic field. In Fig.|5l the dotted line represents the poten- 
tial evolution of the magnetic field, while the dashed line 
represents the dynamical evolution. The important infor- 
mation contained in this diagram can be condensed to the 
following expressions for the two reconnection rates: 
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where Vd cos is the driving velocity perpendicular to the 
direction of the overlying magnetic field, which also rep- 
resents the direction perpendicular to the developing cur- 
rent sheet. This shows that for slow driving velocities 
the dynamic peak reconnection rate approaches the po- 
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Figure 6. Left, the pointing flux for five experiments with 
different driving velocity. The graphs are scaled accord- 
ing to Eq. to give the same absolute value as a func- 
tion of advection distance. Right, the average Joule dis- 
sipation. As a first approximation, the energy released 
from all the experiments, over a given advection distance 
of the sources, may be considered the roughly same. The 
colours of the lines represent different driving velocities. 



tential rate, while for increasing driving velocities it ap- 
proaches a constant rate whereas the potential rate con- 
tinues to grow linearly. For realistic driving velocities in 
the solar context, this indicates that the magnetic recon- 
nection rate, here, is fast, at just a factor of two less than 
the potential rate. 



5. ENERGY CONSIDERATIONS 

Taking the result from the previous section, namely that 
the reconnection is fast and within a factor of two of 
that the required for the field to maintain a potential 
state, one may assume that the energetics of the sys- 
tem could simply be repres ented using a potential model 
jnalsgaardfcParnell 2004). 

If the driving velocity is slow and has a duration that 
is long compared to the Alfven crossing time then the 
Poynting flux input can be estimated by. 
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where Bn is the normal component of the magnetic field, 
Vd is the driving velocity, td is the duration of the driv- 
ing and L is the distance to the anchored (fixed) end of 
the field line. This reveals that the energy input scales 
quadratically with driving speed and linearly with dura- 
tion. Integrating Eq. (|3j over a fixed advection distance 
suggests that the same amount of energy should be in- 
jected to the magnetic configuration independent of the 
driving velocity. These results have been tested by con- 
ducting a series of experiments with different driving ve- 
locities. Scaling the amplitudes of the Poynting flux to 
the advection distance indicates that the results follow the 
simple scalings obtainable from Eq. (|3} to within a veloc- 
ity dependent initial offset, left panel of Fig.|6l The initial 
difference is due to the change in ratio between the driv- 
ing velocity and the local Alfven velocity and holds as 
long as the information can propagate freely away from 



the driving boundary. It is only when the initial perturba- 
tion of the field lines, which is generated when the driving 
is first switched on, reaches another anchor point of the 
field line that the effect of the angle of the field (caused 
by the driving) start to be felt leading to a linear increase 
in time of the Poynting flux according to Eq. (|3}. 

If the same amount of energy is injected into the domain 
for a given advection distance regardless of driving speed, 
then in each case the same amount of energy is available 
for heating the plasma through dissipation. The dynami- 
cal evolution of the various cases follow similar patterns 
and so it is not unreasonable to expect the energy release 
to be of comparable magnitude. The Joule dissipation for 
the numerical domain is plotted as a function of advection 
distance and scaled to the same driving time in right panel 
of Fig.|6l This panel shows that only minor differences in 
the energy release exist between the different cases over 
the same advection distance, and that most of these can be 
accounted for by the differences in the ratio between the 
Alfven crossing time and the driving time for a given ad- 
vection distance. The effect of the energy release, on the 
other hand, depends linearly on the timescale over which 
it is released. The shorter the timescale the greater the 
impact of the energy release. 

The amount of Poynting flux and Joule dissipation can- 
not be estimated using potential models, as they depend 
critically on the magnetic fields reaction to the imposed 
driving, and through this on the change in angle of the lo- 
cal magnetic field lines stressed by the boundary driving. 
This angle is a result of a balanced between the geometry 
of the magnetic field, the local current concentrations and 
the efficiency of the magnetic resistivity. 

6. PLASMA STRUCTURE OF THE ENERGY 
RELEASE 

The models used in these investigations have been setup 
to investigate the effects of magnetic reconnection in a 
simple magnetic structure. But they can also be used to 
attempt to investigate the change in the visual appearance 
of heated magnetic structures as the type of reconnection 
changes with time. The best and only realistic way to do 
this is to include: i) a more detailed model atmosphere; 
ii) the effects of optically thin radiation and anisotropic 
heat conduction in the energy equation. The disadvan- 
tage with this is the requirement for background heating 
to balance the radiation and conduction in the initial con- 
figuration which one cannot easily justify in the given 
system. Instead, a more simplistic approach to illustrate 
the effect has been adopted here. From the field line 
tracing mentioned above it is straight forward to identify 
the starting points of the field lines that change connec- 
tivity between the snapshots. A realistic assumption to 
then make is that the energy released in the reconnec- 
tion event scales with the magnetic flux represented by 
the field line. Distributing this energy evenly along the 
field lines, a local temperature increase of the plasma may 
then be estimated with the inclusion of a vertical den- 
sity profile along the field lines. Here, the plasma density 
is assumed to decreases rapidly with height following a 



Figure 7. The four panels show the structure of the heated 
plasma at different times during the evolution of the flux 
systems in one experiment. Only the regions with temper- 
atures above the initial coronal temperature are shown. 

simple hyperbolic tan function. This simple model em- 
ulates a photosphere at the driving boundary follow by a 
shallow transition region, while most of the domain con- 
tains a coronal plasma density. Then by following the 
reconnected field lines, local temperature increases can 
be calculated which, when all contributions are summed, 
reveal a few interesting conclusions as the system evolves 
in time (Fig.E). 

• The initial connecting of the magnetic field between 
the two sources gives rise to a compact structure out- 
lined by field lines that link the two sources. 

- Note, that a limitation here is the problem with 
identifying the newly reconnected open field 
lines which are not identified. These coronal 
field lines may, however, be so long that the 
energy distributed along them from the recon- 
nection has only limited impact on rising their 
temperature. 

• When the field starts to disconnect the field lines ex- 
tend from the sources to the open coronal field again 
and the energy is distributed along field lines that 
might reach far away from the reconnection site. 

• Furthermore, the disconnecting process is not 
equally active all across the whole separatrix sur- 
face, but is more stochastic in nature leading to 
several fingers of field lines leading away from the 
sources. 

Taking the simplicity of the approach into account it is 
clear that to fully see the thermal, and therefore observa- 
tional appearance of the reconnecting structure, the prob- 
lem has to be solved consistently including the effects of 
optical thin radiation and anisotropic heat conduction. An 



interesting, but difficult task to handle consistently and 
clearly a problem for future investigations. 

7. SUMMARY 

The numerical experiments discussed in this paper inves- 
tigate a simple situation where two initially unconnected 
flux sources are forces into each other by an imposed ve- 
locity profile on the photospheric boundary. From the 
experiments the foUowing main conclusions arises: 
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In terms of observational appearance of the flux interac- 
tion the applied model is very simple and more detailed 
investigations, including a more reaUstic treatment of the 
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